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ABSTRACT 

During  the  1972  AIDJEX  pilot  study,  mast-mounted  current  meters 
were  operated  continuously  for  30  days  in  the  main  camp  to 
measure  currents  at  ten  depths  between  2  m  and  100  m  below  the 
ice  base.  The  currents  were  plotted  in  three  forms:  relative 
speed  and  direction,  absolute  speed  and  direction,  and  progressive 
vector . 


INTRODUCTION 

The  primary  aim  of  the  AIDJEX  oceanographic  program  is  the  determination 
of  water  drag  on  the  base  of  drifting  pack  ice.  During  the  AIDJEX  pilot 
program  in  1972,  ocean  currents  were  monitored  continuously  at  ten  levels 
between  the  surface  and  100  m.  This  report  describes  the  experimental 
techniques  and  data  reduction  and  summarizes  the  results.  An  accompanying 
article  in  this  Bulletin  applies  these  results  to  water  drag. 

A  secondary  objective  of  the  oceanographic  program  is  to  explore 
basic  processes  in  the  upper  layers  of  the  Arctic  Ocean.  The  current  obser¬ 
vations  presented  here  and  those  from  the  1971  pilot  program  have  contributed 
to  a  better  understanding  of  two  phenomena  in  the  upper  layers:  sub-ice 
boundary  currents  and  subsurface  baroclinic  eddies.  The  sub-ice  boundary 
currents  were  seen  most  clearly  in  the  1971  data  and  are  apparently  convection 
cells  driven  by  brine  production  during  freezing.  Subsurface  eddies  were 
most  clearly  observed  in  the  1972  data  presented  here.  The  eddies  are 
10-20  km  in  diameter  and  are  found  between  depths  of  50  and  300  m.  Speeds 
may  reach  40  cm/sec  at  the  core  level,  150  m.  These  eddies  may  play  an 
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important  part  in  the  horizontal  exchange  of  heat,  momentum,  and  salt 
between  the  Arctic  Ocean  and  adjoining  seas.  Details  of  the  eddies  are 
given  in  a  separate  report  [Hunkins,  1974]. 


DATA  ACQUISITION 

Current  meters  were  suspended  on  inverted  masts  with  their  orientation 
rigidly  fixed  to  that  of  the  floe  on  which  the  station  was  situated.  The 
current  meters  (Bendix/Marine  Advisers  Q-9)  were  of  the  Savonius  rotor  type. 
Although  these  meters  have  a  magnetic  compass  for  direction  reference,  the 
compasses  were  locked  by  a  magnet  attached  to  the  exterior  of  the  instrument. 
Current  direction  was  thus  referenced  to  the  floe  itself.  The  azimuth  of 
the  floe  was  monitored  daily  with  celestial  observations. 

The  threshhold  sensitivity  of  the  current  meters  is  conservatively 
estimated  by  the  manufacturer  as  2.6  cm/sec,  although  our  experience  indi¬ 
cates  that  the  threshhold  may  be  considerably  lower.  The  manufacturer's 
calibration  curve  was  used  for  rotor  turns  versus  water  speed.  It  is  a 
general  conclusion  of  most  tests  on  Savonius  rotor  meters  that  this  is 
satisfactory  for  meters  manufactured  with  adequate  quality  control  so  long 
as  the  bearings  are  free.  The  bearings  were  checked  for  clearance  before 
launching  the  instruments  by  observing  rotor  spindown  in  air.  Directional 
accuracy  is  specified  as  ±3°  by  the  manufacturer.  The  current  meters  weigh 
34  pounds  in  air  and  22  pounds  in  water.  Bolted  fittings  attached  the 
masts  to  the  flanges  on  the  top  and  bottom  of  the  instruments.  The  instru¬ 
ments  were  axially  aligned  with  the  mast,  forming  an  extension  of  it,  so 
that  no  supports  interfered  with  flow. 

Three  masts  were  used  to  support  the  ten  current  meters.  Meters 
were  placed  on  the  shallow  mast  at  levels  of  2,  4,  8,  12,  and  20  m  below 
the  ice  base.  All  current  meter  depths  were  referenced  to  the  bottom  of 
the  ice,  nominally  taken  as  2  m  below  sea  level.  (For  depth  below  sea 
level,  add  2  m  to  indicated  depths.)  Ice  thickness  was  measured  at  four 
locations  by  the  CRREL  group  around  the  current  meter  building.  They  found 
ice  thicknesses  of  2.26,  2.04,  2.11,  and  2.08  m  below  sea  level.  Meters 
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were  suspended  at  30,  40,  and  50  m  on  the  intermediate  mast  and  at  70  and 
100  m  on  the  deep  mast. 

The  shallow  mast,  a  prototype  of  the  ones  to  be  used  in  the  1975 
experiment,  was  made  of  one-inch  stainless  steel  tubing  fitted  with  locking 
swivel  attachments  for  alignment  of  individual  meters.  The  intermediate  and 
deep  masts  were  of  lV'  aluminum  pipe  locked  together  with  collars  and  pins. 
Individual  instruments  were  not  aligned  on  the  aluminum  masts,  but  their 
orientations  relative  to  each  other  were  measured.  Proper  alignment  along 
all  masts  was  assured  by  punch  marks  on  the  tubing  or  pipe  applied  in  a 
lathe  bed. 

The  shallow  and  intermediate  masts  were  installed  through  separate 
ice  wells  2  m  apart  in  the  current  meter  building;  the  deep  mast  was  located 
30  m  away  in  the  Lamont  living  quarters.  The  top  of  each  mast  emerged  to 
eye  level  in  the  buildings  and  was  fitted  with  an  azimuth  plate  for  align¬ 
ment.  Surveys  run  with  a  theodolite  tied  the  mast  alignment  to  the  floe 
azimuth. 

Signals,  carried  from  the  sensors  to  the  surface  by  electrical  cables, 
were  converted  to  analog  voltages  by  signal  conditioners  (Bendix/Marine 
Advisers  S-ll) .  Speed  and  direction  at  all  ten  levels  were  displayed  con¬ 
tinuously  on  panel  meters  in  the  current  meter  building  during  the  experiment. 
The  data  were  recorded  on  paper  strip  charts  with  multipoint  servo  recorders 
which  printed  once  per  minute  on  magnetic  tape.  The  digital  data  acquisition 
system  (Hewlett  Packard  Model  2012D)  included  a  reed  scanner,  integrating 
digital  voltmeter,  coupler,  and  incremental  magnetic  tape  recorder.  The 
signal  conditioners  had  a  response  time  of  about  30  seconds  for  the  speed 
channel  and  10  seconds  for  the  direction,  so  that  the  problem  of  aliasing 
high-frequency  fluctuations  with  a  one-minute  sampling  interval  should  be 
negligible.  The  chart  recorders  operated  from  15  March  to  26  April,  but  a 
failure  in  the  tape  recorder  delayed  the  start  of  digital  recording  to 
28  March.  The  digital  data,  essentially  complete  through  the  inclusive 
dates  of  29  March  to  25  April,  are  the  data  described  here. 

The  entire  system,  except  for  the  current  sensors  themselves,  was 
calibrated  at  the  beginning,  the  middle,  and  the  end  of  the  experiment. 
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Calibration  signals  were  introduced  into  the  signal  conditioners  and  monitored 
on  the  chart  recorders  and  digital  voltmeter.  A  precision  decade  resistance 
box  was  used  to  simulate  direction  changes  of  the  vane  and  potentiometer  in 
the  actual  direction  sensor.  A  square-wave  generator  simulated  output  of 
the  magnetic  reed  switch  in  the  actual  speed  sensor.  The  data  were  corrected 
on  the  basis  of  the  calibrations  during  processing  after  the  field  study. 

Since  the  ocean  currents  are  observed  from  a  moving  platform,  a  know¬ 
ledge  of  ice  drift  is  required  for  complete  interpretation  of  the  results. 

Ice  position  was  monitored  at  irregular  intervals  of  roughly  one  hour  with 
a  U.  S.  Navy  satellite  navigation  system.  The  position  data  were  smoothed 
with  a  Kalman  filter  with  a  response  of  50%  at  a  period  of  14  hours.  The 
slope  of  the  Kalman  filter  is  somewhat  steeper  than  that  of  a  cosine  filter. 
Details  of  the  treatment  of  the  navigation  data  are  given  by  Thorndike  [1973]. 

Smoothed  ice  speed  and  direction  are  shown  in  Figure  1.  Since  these 
data  have  been  filtered,  it  is  of  interest  to  know  if  higher-frequency 


Fig.  1.  Smoothed  ice  speed  and  direction,  AIDJEX  main  camp,  1972. 
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motions  of  the  ice  have  been  suppressed.  To  check  this  point,  positions  were 
taken  over  a  period  of  ten  days  with  an  acoustic  bottom  reference  system  for 
comparison  with  the  satellite  navigation.  The  acoustic  system  yields  highly 
precise  position  relative  to  a  bottom  transponder  at  intervals  of  a  few 
minutes.  Combined  data,  incorporating  both  satellite  and  acoustic  results, 
presumably  represent  ice  motion  of  all  periods  in  all  essential  detail. 

These  combined  data  were  compared  with  the  satellite  data  alone.  The 
variance  of  the  difference  between  the  two  time  series  was  1000  m  hr  for 

the  east  component  and  650  for  the  north.  This  indicates  that  the  ice 
velocities  are  reliable  to  better  than  1  cm/sec. 

Inspection  of  the  data  plot  of  Thorndike  reveals  that  much  of  the 
variance  is  associated  with  oscillations  with  a  period  of  about  12  hours. 

This  is  probably  due  to  inertial  oscillations  of  the  ice,  which  are  suppressed 
by  the  smoothing.  For  some  purposes,  12-hour  means  of  ice  velocity  have 
been  used  to  eliminate  these  oscillations. 


DATA  REDUCTION 

The  data  tapes  written  by  the  digital  data  system  are  formatted  BCD 
tapes  that  can  be  read  by  many  computers.  Since  all  the  current  meter  data 
were  to  be  reduced  on  a  single  computer  (IBM  1130) ,  the  tapes  were  rewritten 
in  unformatted  arrays.  This  reduced  the  amount  of  magnetic  tape  used  from 
two  full  2400-foot  reels  to  about  half  a  reel,  with  considerable  savings 
in  computer  access  time.  These  data  were  expressed  in  volts,  the  original 
field  unit.  The  calibration  results  were  then  used  to  convert  the  data  into 
speed  in  cm/sec  and  degrees  from  true  north.  The  final  tape  record  contains 
data  in  21  channels:  10  speed  channels,  10  direction  channels,  and  a  time 
channel . 

Since  the  data  were  recorded  at  one-minute  intervals  for  almost  a 
month,  there  are  more  than  40,000  data  points  per  channel.  This  number  was 
reduced  considerably  by  filtering  the  data  with  an  80-minute  rectangular 
filter  and  sampling  every  hour.  Each  hourly  value  thus  represents  an 
average  of  the  preceding  and  the  succeeding  40  minutes.  These  filtered  data 
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were  recorded  on  disk  files  to  provide  rapid  random  access.  The  one-minute 
data  have  not  been  used  since  they  were  filtered,  but  the  magnetic  tapes 
are  still  available  for  studying  high-frequency  phenomena. 

Ice  velocity  at  one-hour  intervals  was  also  recorded  on  disk.  These 
data  were  then  added  vectorially  to  the  data  for  currents  measured  relative 
to  the  ice  motion.  The  true  currents  thus  obtained  were  also  written  on 
disk.  The  disk  therefore  contains  three  types  of  data,  each  sampled  hourly: 
ice  velocity,  relative  current  velocity  at  ten  depths,  and  true  velocity  at 
ten  depths.  These  are  the  data  sets  that  have  been  used  to  construct  the 
figures  and  tables  in  this  report. 


DATA  PRESENTATION 

Ocean  currents  in  general  form  a  three-dimensional  vector  field  that 
changes  continuously  with  time.  The  current  meters  used  in  AIDJEX  measured 
only  horizontal  currents.  Vertical  currents  are  usually  several  orders  of 
magnitude  less  than  horizontal  currents  in  the  ocean.  The  AIDJEX  data  are 
thus  two-dimensional  vectors  in  time  at  discrete  levels  in  the  ocean.  The 
location  of  the  measuring  site  varies  with  time  as  the  ice  station  drifts. 

The  station  drifted  generally  westward  at  a  rate  of  about  2  km/day.  At  the 
beginning  of  digital  recording  on  29  March  the  station  was  at  75°03'N 
148°43’W;  at  the  end,  on  25  April,  it  was  at  75°06'N  151°32’W.  The  observa¬ 
tions  are  thus  neither  exactly  in  an  Eulerian  nor  exactly  in  a  Lagrangian 
frame  of  reference.  The  first  would  be  from  a  point  fixed  relative  to  the 
ocean  floor,  and  the  second  would  follow  the  path  of  a  water  particle. 

Since  a  time-varying  vector  field  involves  four  dimensions,  it  cannot 
be  simply  depicted  on  a  two-dimensional  page.  Presentations  of  current  data 
are  compromises,  and  some  presentations  are  better  than  others  for  emphasizing 
particular  aspects  of  current  behavior.  From  the  many  possibilities,  two 
types  of  presentation  were  selected  for  use  here.  The  first,  current  speed 
and  direction,  either  relative  or  absolute,  plotted  against  time,  emphasizes 
such  time-varying  current  features  as  tides,  inertial  oscillations,  and 
storm  effects.  The  second  presentation,  progressive  vector  diagrams  representing 
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the  displacement  of  a  water  particle  provided  that  the  currents  are 
unchanging  horizontally,  emphasizes  the  net  displacement  of  a  water  particle 
or  the  mean  current. 

The  quantities  actually  measured  were  current  speed  and  direction 
relative  to  the  ice  at  ten  levels.  The  relative  currents  are  shown  in 
Figures  2-11.  Note  that  the  times  in  these  diagrams  and  in  all  others  in 
this  paper  are  Alaskan  Standard  Time,  where  AST  =  GMT  -  10  hours.  Direction 
is  referenced  to  true  north. 

The  relative  current  plots  provide  some  indication  of  the  operating 
condition  of  the  instruments.  Any  stickiness  of  the  rotor  or  vane,  for 
example,  will  be  evident.  One  such  fault  was  detected,  not  in  the  field, 
but  in  the  plot  (Fig.  3):  a  sluggish  vane  at  4  m  depth,  whose  variation 
is  much  less  than  that  of  the  instruments  above  and  below  it.  A  fault  like 
this  would  be  obscured  by  any  subsequent  manipulation  of  the  data  to  combine 
speed  and  direction  or  to  incorporate  them  into  other  data.  We  show  the 
4  m  data  here  for  completeness,  but  we  do  not  use  them  in  later  interpreta¬ 
tion.  Another  error,  a  misalignment  of  the  current  meters  at  30  and  40  m, 
was  detected  later  by  comparisons  with  auxiliary  data  from  other  current 
meters.  The  directions  were  corrected  during  data  reduction.  All  of  the 
other  current  meters  appear  to  have  produced  satisfactory  results. 

Relative  currents  are  useful  for  drag  calculations  where  the  velocity 
difference  between  ice  and  water  is  important.  For  other  studies  of  ocean 
currents  it  is  generally  useful  to  remove  the  ice  drift.  In  the  extreme 
case  when  the  ice  is  motionless,  the  relative  observations  are  equivalent 
to  true  ocean  currents.  In  the  other  extreme  case  of  moving  ice  but 
motionless  water,  the  relative  currents  reflect  ice  drift  and  the  current 
meter  acts  as  a  pit  log.  In  general,  however,  both  ice  and  water  are  in 
motion,  and  true  ocean  currents  are  obtained  by  adding  the  ice  velocity  to 
the  relative  current  velocity  with  due  regard  for  their  vectorial  nature. 

In  practice,  the  velocities  are  separated  into  their  north  and  east  compo¬ 
nents  and  the  addition  is  performed  separately  on  each  component.  Plots 
of  true  current  speed  and  direction  are  shown  in  Figures  12-21.  Progressive 
vector  diagrams  for  the  true  currents  are  shown  in  Figures  22-31. 
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The  current  measurements  for  this  study  were  made  within  two  ranges 
of  depth,  each  range  characterized  by  a  particular  current  behavior  and 
physical  regime.  The  first  of  these  ranges,  the  planetary  boundary  layer, 
extends  from  the  surface  to  a  depth  of  about  35  m.  (Salinity  and  temperature 
observations  showed  that  the  mixed  layer  also  extended  down  to  about  35  m 
during  the  observational  period,  so  that  these  two  layers  can  be  considered 
here  as  essentially  the  same.)  In  the  planetary  boundary  layer,  currents 
move  primarily  under  the  influence  of  friction,  horizontal  pressure  gradients, 
and  the  earth's  rotation.  There  is  a  large  shear  in  the  upper  part  (5-10  m) 
of  this  layer  during  storms,  but  little  shear  occurs  at  any  time  in  the 
lower  part,  where  the  flow  is  apparently  nearly  geostrophic.  Since  the 
water  is  also  nearly  homogeneous,  it  is  essentially  a  barotropic  flow. 

The  deep  current  measurements  (at  40,  50,  70,  and  100  m)  were  taken 
in  the  steep  density  gradient  that  extends  downward  from  the  planetary 
boundary  layer  to  a  depth  of  300  m.  In  this  region  the  currents  also  appear 
to  be  in  nearly  geostrophic  equilibrium,  but  the  flow  is  highly  baroclinic, 
with  strong  vertical  shears  at  all  times.  The  swiftest  currents  appear 
at  the  deepest  levels  and  are  associated  with  baroclinic  subsurface  eddies. 
These  eddies  appear  intermittently  on  the  records,  attaining  speeds  as 
great  as  35  cm/sec  at  100  m. 

Means  over  periods  longer  than  one  hour  are  useful  if  currents  are 
reasonably  steady  over  the  averaging  interval  chosen.  The  longer  intervals 
tend  to  smooth  out  effects  of  turbulence  and  oscillating  currents.  A  12-hour 
interval  was  chosen  to  eliminate  the  effects  of  inertial  and  semidiurnal 
tidal  motions,  which  both  have  periods  close  to  12  hours  and  which  are 
noticeable  in  the  ice  velocity  records  and  in  some  of  the  water  velocity 
records.  Both  relative  and  absolute  current  means  over  12-hour  intervals 
are  given  at  all  ten  depths  in  the  Appendix.  Mean  ice  velocity  is  also 
given. 
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Figs.  2-11.  Relative  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  2-11.  Relative  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  2-11.  Relative  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  2-11.  Relative  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  2  11.  Relative  current  speed  and  direction  for  specified  depths  below  the  base  of  the  i 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  12-21.  True  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  12-21.  True  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  12-21.  True  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Figs.  12-21.  True  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Fig.  24.  Depth  =  8  meters. 
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Fig.  25.  Depth  =  12  meters. 
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Fig.  26.  Depth  =  20  meters. 
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Figs.  22-31.  Progressive  vector  diagrams  for  currents  at  specified 
depths  below  the  base  of  the  ice  (sea  level  =  2  m) ,  plotted  at 
hourly  intervals . 
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Figs.  12-21.  True  current  speed  and  direction  for  specified  depths  below  the  base  of  the  ice 
(sea  level  =  2  m) ,  plotted  at  hourly  intervals. 
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Fig.  22.  Depth  =  2  meters 


Fig.  23.  Depth  =  4  meters. 


Figs.  22-31.  Progressive  vector  diagrams  for  currents  at  specified 
depths  below  the  base  of  the  ice  (sea  level  =  2  m) ,  plotted  at 
hourly  intervals . 
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Fig.  27.  Depth  =  30  meters. 
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Fig.  29.  Depth  =  50  meters. 
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Figs.  22-31.  Progressive  vector  diagrams  for  currents  at  specified 
depths  below  the  base  of  the  ice  (sea  level  =  2  m) ,  plotted  at 
hourly  intervals. 
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Figs.  22-31.  Progressive  vector  diagrams  for  currents  at  specified 
depths  below  the  base  of  the  ice  (sea  level  =  2  m) ,  plotted  at 
hourly  intervals. 
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APPENDIX 


CURRENT  OBSERVATIONS  AT  THE  MAIN  CAMP 
1972  AIDJEX  PILOT  STUDY 


12-hour  mean  values 
Velocity  components  in  cm/sec 
Time  is  Alaskan  Standard  Time 
(AST  =  GMT  -  10  hr). 

Direction  is  given  in  terms  of  true  north. 
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THE  OCEANIC  BOUNDARY  LAYER  AND  ICE-WATER  STRESS 
DURING  AIDJEX  1972 


by 

Kenneth  Hunkins 

Lamont- Doherty  Geological  Observatory 
Columbia  University 3  Palisades  3  New  York  10964 


ABSTRACT 

During  the  periods  of  storm  and  rapid  ice  drift,  currents 
in  the  upper  15  to  25  m  showed  development  of  a  modified  Ekman 
spiral  during  the  1972  AIDJEX  pilot  study.  A  momentum  integral 
method  is  used  to  evaluate  ice-water  stress.  Maximum  hourly 
mean  stress  is  1.66  dyn/cm2.  The  balance  of  forces  on  the  ice 
for  cases  when  wind  speed  exceeds  5  m/sec  shows  a  consistent 
pattern.  Ice-water  stress,  air-water  stress,  and  Coriolis  force 
are  all  of  the  same  magnitude.  The  pressure  gradient  force  is, 
however,  much  smaller  than  these.  Internal  ice  resistance  is 
found  as  a  residual  from  the  force  diagram  and,  during  the 
rapid  drift  periods,  is  directed  about  135°  to  the  left  of  the 
ice  drift. 


The  exchange  of  properties  between  the  polar  air  and  the  polar 
water  takes  place  across  the  pack  ice  barrier  moving  on  the  interface 
between  atmosphere  and  ocean.  In  the  case  of  momentum  exchange,  the  flow 
is  generally  downward  from  air  to  water.  Wind  stress  on  the  upper  ice 
surface  acts  as  a  driving  force  and  water  stress  on  the  lower  ice  surface 
acts  as  a  retarding  force.  Momentum  is  transported  downward  from  the 
atmosphere  through  the  ice  and  into  the  oceanic  planetary  boundary  layer. 

One  of  the  aims  of  the  oceanographic  program  during  the  AIDJEX  pilot 
experiment  in  1972  was  to  assess  this  momentum  exchange  and  measure  ice-water 
stress  in  the  central  Arctic  Ocean.  For  this  purpose  detailed  current 
observations  were  made  in  the  upper  levels  of  the  ocean  from  the  lower  ice 
surface  through  the  planetary  boundary  layer  and  below  it.  The  structure 
of  the  boundary  layer  and  its  changes  with  time  over  a  period  of  one  month 
were  recorded.  From  this  information  the  ice-water  stress  was  calculated. 
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The  techniques  of  current  measurement  and  data  reduction  are  described 
in  an  accompanying  article  in  this  issue  of  the  AIDJEX  Bulletin.  Data 
averages  over  intervals  of  12  hours  are  used  here  to  minimize  fluctuations 
due  to  turbulence  and  inertial  oscillations.  Over  intervals  of  12  hours 
ice  drift  and  ocean  currents  are  generally  quite  stable,  so  that  we 
feel  justified  in  treating  these  as  equilibrium  situations.  Winds 
were  light  and  ice  drift  small  during  much  of  the  four  weeks.  Most  of  the 
drift  was  driven  by  the  high  winds  of  short-lived  storms.  Only  selected 
intervals  when  ice  and  wind  motion  was  sufficiently  high  were  used  for  study 
here.  During  low  winds  the  pattern  is  not  so  evident  and  is  obscured  by 
random  experimental  errors,  but  during  high  winds  clear  patterns  of  current 
structure  emerge. 

Current  hodographs,  lines  connecting  the  tips  of  current  vectors,  are 
presented  in  Figures  1  and  2  for  intervals  when  mean  ice  drift  exceeded 

V,  CM/SEC 


Fig.  1.  Hodographs  of  absolute  currents  for  12-hour  intervals  with  ice  drift 
greater  than  9  cm/sec.  Depths  in  meters  below  ice  base.  North  at  top  of 
page.  Local  standard  time. 
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9  cm/sec.  A  mean  hodograph  over  a  60-hour  period  coinciding  with  the 
largest  storm  is  shown  in  Figure  3.  The  zero  depth  level  for  these  measure¬ 
ments  is  the  base  of  the  ice,  2  m  below  sea  level.  The  hodograph  is  dashed 
in  the  upper  2  m;  this  is  considered  to  be  the  surface  boundary  layer  in 
which  frictional  forces,  but  not  the  earth's  rotation,  are  important.  In 
the  surface  boundary  layer  there  is  a  large  vertical  current  shear  in  speed 
with  very  little  change  in  current  direction.  Below  the  2  m  level  there  is 
a  large  current  shear  in  both  speed  and  direction.  This  is  the  Ekman  layer, 
in  which  the  effect  of  the  earth’s  rotation  is  manifested  by  a  clockwise 
rotation  of  the  current  direction.  A  theoretical  Ekman  spiral  decreases 
exponentially  downward  in  speed,  the  angle  increasing  in  direct  proportion 
to  depth.  The  spiral  in  some  of  the  plots  tends  to  continue  downward  to 
50  m.  At  the  center  of  the  decreasing  spiral  is  the  geostrophic  current 
vector  at  the  base  of  the  Ekman  layer  where  frictional  effects  vanish. 
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Fig.  2.  Hodographs  of  absolute  currents  for  12-hour  intervals  with  ice 
drift  greater  than  9  cm/sec.  Depths  in  meters  below  ice  base. 
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Fig.  3.  Mean  current  hodograph  for  60-hour  interval  during  storm. 

The  exact  depth  of  an  Ekman  spiral  is  not  well  defined  since  it  spirals 
infinitely  downward.  Ekman  took  the  depth  to  be  that  at  which  the  current 
flowed  in  the  direction  opposite  to  that  at  the  surface.  Others  have  chosen 
the  level  at  which  the  current  speed  has  decreased  by  a  factor  of  1/e  from 
its  surface  value.  The  exact  choice  is  not  critical  for  most  purposes, 
since  the  speed  decreases  exponentially  and  changes  are  small  near  the  base 
of  the  layer.  The  geostrophic  current  vector  would  lie  in  the  center  of  the 
loop  between  25  and  50  m. 

The  mixed  layer  revealed  by  temperature  and  salinity  profiles  during 
the  experiment  was  about  35  m  deep.  Below  this  level  a  steep  gradient  in 
these  properties  and  density  extended  to  a  depth  of  300  m.  The  upper  part 
of  the  mixed  layer,  0-15  m,  was  often  unstable,  with  brine-driven  convection 
caused  by  freezing  and  brine  production  at  the  ice  base  [Smith,  in  press]. 

The  lower  part  of  the  mixed  layer,  15-35  m,  was  generally  neutral  or  slightly 
stable.  The  Ekman  layer  and  frictional  exchange  of  momentum  must  be  limited 
by  the  sharp  density  interface  below  35  m  and  to  some  extent  by  the  change 
at  about  15  m. 
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MOMENTUM  INTEGRAL  METHOD 


A  technique  involving  direct  evaluation  of  momentum  in  the  upper  water 
column  from  current  profiles  offers  several  advantages  for  determining  ice- 
water  stress  on  pack  ice.  Effects  of  both  friction  and  form  drag  are 
included  in  the  result.  The  assumptions  are  few,  and  none  is  made  about 
the  nature  of  eddy  viscosity.  The  method  can  be  demonstrated  most  simply 
for  the  case  of  smooth  ice  so  that  the  nonlinear  and  pressure-gradient  terms 
are  neglected,  i.e.,  only  skin  friction  is  acting.  In  the  appendix  it  is 
shown  that  the  method  is  more  general  and  that  with  a  proper  set  of  current 
observations  the  form  drag  can  also  be  included  in  the  stress  result.  The 
horizontal  equations  of  motion  in  their  vertically  integrated  form  are 
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where  the  mass  transports  are  given  by 

f  0 

Mx  =  J  P(w  -  u  )dz 

-H 

My  =  /  P(y  "  vg)dz 

~H 


with  z  positive  upwards.  It  is  assumed  that  there  is  a  stress  on  the 

./  V 

planetary  boundary  layer  with  components  T  and  T  at  the  ice-water 

o  tx  o  ,zy 

interface  but  no  stress  at  the  base,  of  the  layer.  Observed  velocity 

components  are  uy  V,  while  w  ,  Vg  are  the  geostrophic  velocity  components 
at  depth  -H,  where  frictional  effects  become  negligible. 

The  Cartesian  system  used  here  is  a  valid  approximation  for  small 
areas  at  some  distance  from  the  pole.  Very  close  to  the  pole  a  spherical 
coordinate  system  would  be  necessary.  A  constant  Coriolis  parameter,  /, 
is  also  a  valid  approximation  for  the  relatively  small  areas  considered. 
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difference.  For  comparison,  wind  stress  on  the  upper  ice  surface  was 
calculated  from  the  synoptic  observations.  The  drag  law 


was  used  with  p  =  0.00125  g/cm3  and  C  =  1.5  x  10  3.  The  drag  coefficient 

Si  L) 

is  based  on  eddy  correlation  measurements  made  with  a  sonic  anemometer  by 
the  Bedford  Institute  of  Oceanography  group.  Because  of  the  quadratic 
nature  of  the  drag  law  it  is  not  strictly  appropriate  to  use  the  average 
wind  speeds.  In  this  case,  however,  their  use  seems  justified  since  the 
observed  winds  were  fairly  stationary  over  the  12-hour  time  intervals.  The 
reverse  of  the  wind  stress  direction  is  plotted  in  Figure  4  for  easier 
comparison  with  water  stress  direction;  the  stresses  are  generally  nearly 
opposed  to  each  other. 

Peak  stresses  on  both  top  and  bottom  of  the  ice  occurred  during  the 
storm  of  12-13  April.  It  is  interesting  that  there  is  a  tendency  for  the 
water  and  wind  stress  to  track  each  other  fairly  closely  for  the  first  day 
or  two  of  a  storm.  After  that,  water  stress  decreases  although  wind  stress 
continues  to  climb  for  another  12  hours.  This  pattern  is  also  evident  in 
the  storm  of  17-18  April.  Maximum  mean  wind  stress  over  12  hours  exceeded 
the  water  stress,  reaching  a  maximum  value  of  1.875  dyn/cm2 ,  which  coincided 
with  a  wind  speed  of  10.0  m/sec  between  1200  and  2400  LST  on  14  April.  The 
water  stress  maximum  of  1.101  dyn/cm2  was  attained  in  the  preceding  12-hour 
interval. 


BALANCE  OF  FORCES 

During  periods  of  strong  winds,  swift  ice  motion,  and  high  stresses, 
conditions  are  frequently  fairly  stationary  over  12-hour  intervals  and  the 
forces  acting  may  be  considered  to  be  in  equilibrium.  Information  on  the 
surface  stresses  of  wind  and  water  has  already  been  presented.  The  body 
forces  due  to  the  pressure  gradient  associated  with  tilt  of  the  ocean 
surface  and  to  the  Coriolis  force  associated  with  the  earth’s  rotation  must 
be  determined.  Since  the  current  at  25  m  has  been  assumed  to  be  in  geos trophic 
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Fig.  4.  Ice-water  (solid  line)  and  air-ice  (dashed  line)  stress  at  1972 
AIDJEX  main  camp  based  on  12-hour  means.  Direction  of  air-ice  stress 
reversed  for  comparison. 
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The  variation  of  this  parameter  with  latitude  would  have  to  be  included  if 
large  areas  of  the  earth’s  surface  were  being  considered. 

The  geostrophic  current  is  assumed  to  be  barotropic,  since  the  planetary 
boundary  layer  coincided  with  the  upper  mixed  layer  in  the  1972  experiments 
so  that  no  baroclinic  currents  would  have  been  present.  A  uniform  geostrophic 
current  prevails  through  the  Ekman  layer.  Its  value  is  the  observed  current 
at  depth  -H .  At  deeper  levels  there  is  stratification  and  there  may  be 
vertical  shear  in  the  geostrophic  current.  At  shallower  depths  the  observed 
current  is  the  sum  of  frictional  and  geostrophic  currents.  If  density  varia¬ 
tions  are  present  in  the  upper  layer,  hydrographic  stations  would  be  necessary 
to  evaluate  the  baroclinic  contribution  to  the  geostrophic  current,  but  for 
a  well-mixed  Ekman  layer,  current  observations  alone  are  sufficient  to  deter¬ 
mine  ice-water  stress  by  the  momentum  integral  technique. 

The  method  depends,  like  any  experimental  method,  on  the  quality  of 
the  observations.  Sufficient  measurements  must  be  made  of  currents  in  a 
vertical  profile  and  spaced  closely  enough  in  time  that  representative  means 
can  be  formed.  The  available  data  described  in  the  accompanying  article 
appear  to  be  adequate  for  the  purpose.  From  the  data  a  choice  of  depth,  -H, 
must  be  made.  It  is  evident  in  Figures  1  and  2  that  a  mean  Ekman  spiral 
exists  over  time  periods  of  12  hours.  The  depth  -H  may  be  arbitrarily 
defined  as  the  depth  at  which  the  current  flows  opposite  to  its  surface 
direction,  or  it  may  be  taken  as  the  depth  at  which  speed  has  decreased  to 
1/e  of  its  surface  value.  Another  approach  would  be  to  choose  the  depth  as 
the  first  minimum  in  the  north  and  east  velocity  profiles.  At  that  level, 
mean  stress  would  vanish.  Still  another  approach  is  to  examine  the  hydro- 
graphic  profiles  for  the  sharp  density  change  at  the  base  of  the  mixed  layer 

and  assume  that  this  coincides  with  the  depth  of  the  Ekman  layer.  Any  one 

or  some  combination  of  these  ways  of  choosing  -H  could  be  objectively  stated 

so  that  the  level  could  be  selected  individually  for  each  profile.  This  may 

be  done  in  the  future,  but  for  the  1972  data  a  constant  level  of  25  m  was 
chosen  for  all  profiles.  This  is  approximately  the  level  of  current  reversal 
and  decrease  of  speed  to  1/e  of  the  value  at  2  m  relative  to  the  geostrophic 
current.  Since  speeds  are  small  near  this  level,  the  exact  choice  of  depth 
does  not  seem  especially  critical. 
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An  added  guide  to  the  choice  of  Ekman  layer  depth  is  provided  by 
hydrographic  profiles  which  show  instability  due  to  brine  convection  in  the 
0-15  m  layer,  neutral  or  slightly  positive  stability  in  the  15-35  m  level, 
and  strong  stability  below  35  m.  Convection  in  the  uppermost  15  m  will  mix 
momentum  thoroughly  down  to  this  depth  and  possibly  somewhat  below  it.  The 
choice  of  25  m  here  is  in  agreement  with  these  hydrographic  measurements. 

The  equations  imply  horizontal  uniformity,  and  only  a  single  vertical 
current  meter  array  was  used.  Keels  of  pressure  ridges  may  be  expected  to 
extend  to  an  average  depth  of  10  m  at  an  average  spacing  of  about  200  m. 

Averaging  over  these  keels  is  accomplished  by  the  horizontal  transport  of 
momentum  in  the  Ekman  layer.  Effects  of  keels  upstream  are  carried  to  the 
observation  site  by  advection.  The  distance  of  influence  may  be  loosely 
estimated  as  one  to  two  orders  of  magnitude  of  the  Ekman  layer  thickness, 
that  is,  250  m  to  2500  m.  Another  estimate  of  distance  of  upstream  influence 
can  be  found  from  the  response  time  of  an  Ekman  layer,  which  is  somewhat 
shorter  than  one  inertial  period  (12.4  hr  at  latitude  75°N).  At  10  cm/sec 
a  water  particle  travels  a  distance  of  about  2  km  in  one  inertial  period. 

Ice-water  stress  was  evaluated  for  the  period  29  March  to  25  April 
using  hourly  current  values  at  levels  of  2,  8,  12,  20,  and  30  m  below  the 
base  of  the  ice.  Differences  of  current  velocity  between  the  geostrophic 
level  at  25  m  and  shallower  levels  were  calculated  and  integrated  vertically 
using  the  trapezoidal  rule.  These  are  the  Ekman  mass  transports;  they 
assume  a  water  density  of  1.0  g/cm3 .  The  time-dependent  term  was  found  by  taking 
hourly  differences  of  these  transports.  The  mean  time-dependent  term  was 
about  one  order  of  magnitude  less  than  the  Coriolis  force  term.  Note  that 
only  relative  currents  are  required  for  finding  Ekman  mass  transport. 

Variations  of  stress  magnitude  and  direction  are  shown  in  Figure  4. 

The  maximum  mean  stress  over  12  hours  was  1.10  dyn/cm  for  0000-1200  LST 
on  12  April.  The  maximum  mean  hourly  stress  was  1.66  dyn/cm2  at  0800-0900  LST 
of  the  same  day.  Stress  results  during  the  same  period  were  reported  as 
1.0  dyn/cm  by  the  University  of  Washington  using  an  eddy  correlation  method 
[Smith,  in  press].  The  eddy  correlation  method  gives  only  surface  friction, 
and  it  may  be  that  the  difference  between  the  two  methods  is  attributable 
to  form  drag.  However,  the  errors  in  the  methods  may  account  for  the 
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balance  and  since  the  pressure  gradient  does  not  vary  with  depth  over  the 

top  25  m  of  the  water  column,  the  pressure  gradient  found  at  the  base  of 

the  Ekman  layer  is  the  same  as  that  acting  on  the  ice.  F  „  =  -p.  h f V „  and 

p,E  1  y 

F  =  p  h f u„  are  the  pressure-gradient  forces  in  the  east  and  north  direc- 

p,N  i  y 

tions,  respectively,  where  the  following  constants  are  used:  ice  density, 
p.  =  0.9  g/cm3;  ice  thickness,  h  =  2.5  m;  and  Coriolis  parameter,  f  = 

2ft  sin4>,  where  cj)  =  75°  and  is  the  earth's  angular  velocity. 

Note  that  absolute  current  velocities  are  required  to  find  the  pressure- 
gradient  force.  The  Coriolis  force  depends  only  on  ice  velocity  and  is 
given  by 


P±hfV 


P  ±hfU 


where  U  and  V  are  components  of  ice  velocity  in  the  east  and  north  directions, 
respectively . 

These  four  forces  are  plotted  for  selected  12-hour  intervals  in 
Figures  5,  6,  and  7  and  for  60-hour  intervals  in  Figure  8.  There  is  also 
an  inertial  force  due  to  the  acceleration  of  the  ice, but  this  term  is 
negligible  in  the  present  cases,  being  too  small  to  plot.  The  eight  cases 
in  these  figures  are  for  intervals  when  mean  wind  speed  exceeded  5  m/sec. 

With  the  constants  used,  this  is  equivalent  to  cases  when  wind  stress 
exceeded  0.5  dyn/cm  .  The  velocity  of  the  ice  is  plotted  also  in  these 
diagrams,  although  it  does  not  enter  into  the  balance  of  forces.  If  it  is 
assumed  that  the  ice  is  drifting  under  balanced  forces,  then  an  equilibrant 
is  required  for  balance.  This  fifth  force,  not  measured  directly  but 
determined  as  a  residual,  is  the  internal  ice  force. 

The  diagrams  indicate  that  all  forces  except  the  pressure-gradient 
are  generally  significant.  However,  on  longer  time  scales,  pressure-gradient 
forces  may  assume  more  significance.  Over  a  mean  drift  of  one  month  for 
AIDJEX  1972  the  Coriolis  and  pressure-gradient  forces  are  nearly  balanced 
[Newton  and  Coachman,  1973] . 
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Fig.  5.  Balance  of  forces  on  ice  for  12-hour  intervals  when  winds  exceeded  5  m/sec. 

Ice  velocity,  V.  ,  in  cm/sec;  air-ice  stress,  T^;  ice-water  stress,  Ty;  Coriolis  force 
F  ;  pressure-gradient  force,  F  ;  and  internal  ice  force,  ^int* 
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Fig.  6.  Balance  of  forces  on  ice  for  12-hour  intervals  when  winds 
exceeded  5  m/sec.  Ice  velocity,  in  cm/sec;  air-ice  stress,  Ta; 

ice-water  stress,  tw;  Coriolis  force,  Fc;  pressure-gradient  force, 
Fpj  and  internal  ice  force,  F^nt_. 
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Fig.  7.  Balance  of  forces  for  two  cases  with  westerly  winds. 
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Fig.  8.  Balance  of  forces  for  60-hour  interval. 


The  pattern  of  stresses  is  consistent  during  the  highest  values  of 
stress  (Fig.  5).  During  some  lower  values  of  stress  (Fig.  6),  deviations 
from  the  earlier  pattern  are  seen.  The  consistency  of  the  pattern  during 
the  highest  stresses  suggests  that  random  errors  are  not  a  problem  and  that 
averaging  in  time  and  space  is  adequate.  The  presence  of  systematic  errors 
is  not  ruled  out,  however.  All  of  these  were  cases  of  easterly  winds.  Two 
cases  of  westerly  winds  were  examined  to  see  if  the  pattern  changes  with 
azimuth.  These  cases  (Fig.  7)  were  for  winds  less  than  5  m/sec.  The  force 
diagram  still  shows  a  similar  pattern  for  these  two  cases,  with  air  stress 
about  30°  to  the  left  of  ice  motion,  water  stress  about  135°  to  the  right 
of  ice  motion,  and  internal  ice  stress  about  135°  to  the  left  of  ice  motion. 
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APPENDIX 

THE  INTEGRATED  STRESS  ON  A  ROUGH  SURFACE  IN  A  ROTATING  SYSTEM 


Consider  a  rough  surface,  such  as  the  underside  of  pack  ice,  moving 
over  an  ocean  of  great  depth.  Frictional  effects  will  vanish  at  H,  the 
depth  of  the  planetary  boundary  layer,  which  in  this  case  is  small  relative 
to  the  total  depth  of  the  ocean.  The  ^-component  of  the  Navier-Stokes 
equation  in  a  local  Cartesian  coordinate  system  is 
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where  u,  V,  and  W  are  the  velocity  components  in  the  x,  y,  and  s  directions, 

respectively;  p  is  density;  p  is  pressure;  and  T  ,  T  ,  and  T  are  the 

CCCC  £jCC 

viscous  stress  components.  The  equation  of  continuity  is 


3 u  3y  3w 
3x  3  y  3s 


(2) 


for  uniform  density  conditions.  We  are  concerned  with  flows  on  scales  small 
enough  that  we  may  neglect  the  variation  of  /,  the  Coriolis  parameter,  with 
latitude.  We  are  interested  in  average  stress  over  a  distance,  X,  sufficiently 
long  to  include  a  representative  sample  of  roughness  conditions.  A  weighting 
factor  is  defined  to  average  horizontally  across  the  water  between  solid 
protuberances.  Average  values  in  the  ^-direction  of  a  representative  depend¬ 
ent  variable,  q,  are  defined  as 

XI 2 

q(y,z9t)  =7  /  a  q (x)dx  (3) 

-X/2 

The  weighting  factor  a  is  unity  when  the  position  x  lies  within  the  fluid 
and  is  zero  when  it  lies  within  a  solid  portion  of  the  rough  boundary. 
Perturbations  from  the  horizontal  average  are  defined  by 

q'  =  q  -  q  (4) 

Substituting  equation  3  into  equation  1  gives  the  perturbation  equation 
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This  equation  is  now  averaged  over  x,  noting  that  the  ^-average  of  an 
^-derivative  gives  values  at  the  points  x  ^  and  x ^  (Fig.  9).  For  example, 
the  average  of  3p'/3x  is  given  by 
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For  z  <  hy  the  greatest  protuberance  depth, 
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the  partial  x  distance, 

), 


is  the  distance  occupied  by  water. 


z=o 
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z  =  -H - 


Fig.  9.  Scheme  for  horizontal  integration  at  constant 
height,  z,  passing  through  solid  keels.  The  weighting 
factor  a  is  zero  within  the  solid  and  unity  in  the 
water. 
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The  x-average  of  the  nonlinear  terms  in  eq.  4  can  be  rewritten  as  in 
the  following  example: 


(V  +  v')—(u  +  u’) 


-du  ,  du  -  du  ’  ,  du  ’ 

dy  ay  dy 

-  du  .  —  du  -  du’  ,  du' 

=  y  -r—  +  V  '  3—  +  V  - -  +  V  '  - - 

dy  dy  dy  dy 

-  du  .  du ’ 

=  V  —  +  V  '  - - 

dy  dy 


(7) 


since  u'  =  v'  =  0  by  definition.  Thus  the  left-hand  side  of  equation  4, 
averaged  over  x,  becomes 
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Substitution  of  the  perturbation  velocity  values  into  the  equation  of 
continuity  (2)  shows  that  the  last  term  in  parentheses  in  eq.  8  is  zero. 

Equation  5  in  ^-averaged  form  then  becomes 
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The  large-scale  pressure  gradient  has  been  separated  out  as  the  geostrophic 
velocity,  u  .  The  remaining  pressure  term  is  associated  only  with  form  drag 

Note  that  the  no-slip  boundary  condition  requires  that  the  term 

h  Cu'u,)  =  °- 
1  ^ 

The  —  Ap '  term  represents  the  net  pressure  drag  across  the  protuber- 
P 

ances,  and  the  Ax^x  term  represents  the  net  effect  of  asymmetries  in  the 
fluctuating  gradient  du’/dx  evaluated  at  positions  x^n  and  cc^.  The  last 
two  terms  on  the  right-hand  side  of  eq.  9  are  the  horizontally  integrated 
effects  of  skin  friction,  including  both  viscous  and  turbulent  contributions 
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Assume  now  that  the  variation  in  the  ^-direction  of  the  average 

parameters,  u,  T  ,  and  v'u\ is  negligible.  This  implies  that  conditions 

yx 

are  laterally  uniform.  Also  assume  that  W  =  0.  Vertical  motion  is  zero 
at  the  upper  boundary  surface  which  is  rigid,  and  this  condition  implies 
that  vertical  motion  vanishes  also  at  the  base  of  the  planetary  boundary 
layer.  Ekman  divergence  effects  which  can  induce  vertical  velocity  at  the 
base  of  the  layer  are  thus  neglected. 

Equation  9  is  integrated  vertically  from  the  base  of  the  planetary 
boundary  layer,  -H ,  where  geostrophic  velocities,  Ug  and  V  ,  exist  upward 
to  the  rough  surface: 
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where  X  is  the  horizontal  weighting  factor  and  the  horizontal  mass  trans¬ 
ports  are  given  by 
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X  (u  -  un)dz 
n  y 
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X  (v  -  vn)dz 
n  9 


Only  the  velocity  departure  from  the  geostrophic  is  used  in  the  case  of 
A/  as  well  as  M .  Since  the  time  rate  of  change  of  M  is  involved,  adding 
a  nearly  constant  geostrophic  velocity  does  not  affect  the  result.  Since 
density  is  constant,  the  geostrophic  velocities  do  not  vary  with  depth. 

The  assumption  of  vanishing  stress  at  the  depth  -H  allows  the  lower 

limit  of  integration  in  the  second  integral  to  be  reduced  from  -H  to  -h. 

The  right-hand  side  of  eq.  10  sums  the  contributions  from  turbulent  and 

roughness  asymmetries,  pressure  drag,  and  skin  friction  due  to  both  molecular 

viscosity  and  turbulent  contributions.  These  effects  may  be  combined  into 

a  single  integrated  boundary  stress  term  with  components  T  and  T 

o  ,x  o 

so  that  eq.  (11)  becomes 
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(12) 
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where  the  z/-equation  has  been  obtained  in  an  analogous  manner  to  the 
x-equation.  Thus  the  total  stress  on  a  rough  surface  in  a  rotating  system 
may  be  determined  from  the  mass  transports  in  the  planetary  boundary  layer 
which  have  been  averaged  over  a  representative  horizontal  distance. 

Equation  12  is  a  statement  of  the  Ekman  transport  relation  including  time 
dependence  and  with  the  added  generalizaiton  that  the  stress  term  includes 
both  skin  friction  and  form  drag.  Vertical  integration  has  freed  the  result 
from  dependence  on  the  details  of  flow  structure  and  friction  within  the 
Ekman  layer. 


This  derivation  is  similar  to  one  given  by  Faller  and  Mooney  [1971] 
for  Ekman  boundary  layer  stress  over  a  rough  surface  in  a  rotating  tank. 


Application  of  the  method  requires  that  the  mass  transports  be 
averaged  over  a  representative  area.  Provided  that  the  averaging  is 
adequate,  the  method  will  be  rigorous  even  though  the  protuberances 
penetrate  a  large  fraction  of  the  planetary  boundary  layer.  In  practice, 
measurements  at  only  a  single  location  have  been  used  for  measuring  stress 
below  pack  ice.  This  is  justified  provided  that  the  site  is  located  at  an 
average  ice  thickness  so  that  no  weighting  of  the  velocity  values  is  required 
and  provided  that  the  protuberances  are  not  so  deep  as  to  cause  extreme 
variations  from  one  location  to  another.  Camp  sites  are  normally  chosen 
away  from  large  ridges  in  areas  of  fairly  level  ice;  this  choice  is  suitable 
for  the  method.  Details  around  large  ridges  are  ignored,  and  for  this  reason 
the  momentum  integral  method  tends  to  underestimate  drag. 
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ABSTRACT 

As  pack  ice  drifts  in  a  stratified  ocean,  the  deep  keels 
of  pressure  ridges  generate  internal  waves  which  produce  a 
retarding  force  on  the  motion.  An  estimate  of  this  drag  force 
was  made  by  scaling  results  from  experiments  performed  long  ago  by 
Ekman,  who  towed  ship  models  in  a  tank  with  two  fluid  layers  of 
different  densities.  For  conditions  assumed  to  be  typical  of 
the  Arctic  Ocean,  internal  wave  drag  is  shown  to  amount  to  only 
10%  of  the  form  drag  and  only  20%  of  the  skin  friction.  The 
two-layered  model  gives  only  the  effect  of  one  mode.  Further 
studies  with  multiple  layers  or  with  continuous  stratification 
would  be  valuable  to  find  the  drag  contribution  of  higher  modes. 


The  resistance  of  fluids  to  moving  bodies  can  usually  be  divided  for 
descriptive  and  measurement  purposes  into  three  types.  Skin  drag  is  caused 
by  the  velocity  shear  near  the  surface  of  the  body.  Form  drag  is  due  to 
the  turbulent  wake  which  occurs  after  the  boundary  layer  has  separated  from 
the  surface.  Both  of  these  types  of  drag  are  a  function  of  velocity,  usually 
squared,  with  a  proportionality  constant  known  as  the  drag  coefficient.  Both 
types  of  drag  have  been  discussed  and  measured  to  some  extent  on  pack  ice. 

If  the  fluid  has  a  free  surface,  a  third  type  of  resistance,  wave  drag, 
may  also  exist.  The  loss  of  energy  to  the  train  of  gravity  waves  left  behind 
the  object  produces  drag.  Surface  wave  drag  is  not  important  in  pack  ice 
where  nearly  total  ice  cover  prevents  the  generation  of  waves  on  the  surface. 
However,  if  a  fluid  is  composed  of  layers  with  differing  densities,  internal 
waves  may  be  generated  that  can  also  produce  wave  drag.  Like  most  oceans, 
the  Arctic  Ocean  does  possess  density  stratification.  Since  internal  waves 
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and  pack  ice  both  have  similar  speeds  it  is  possible  that  the  movement  of 
pack  ice  with  deep  keels  might  generate  significant  internal  waves.  Little 
attention  has  been  given  to  this  type  of  drag  on  pack  ice.  This  note  attempts 
to  provide  an  estimate  for  wave  drag  on  pack  ice  by  scaling  the  results  of 
model  experiments . 

The  first  studies  of  internal  wave  drag  were  prompted  by  the  "dead 
water"  resistance  met  by  ships  under  certain  conditions.  The  effect  of 
internal  wave  drag  on  ships  is  evident  at  low  speeds  of  about  one  knot,  and 
in  areas  where  a  thin  layer  of  fresh  water  overlies  salt  water.  Thus  it  was 
especially  noticeable  in  the  time  of  slow  sailing  ships  and  in  areas,  such 
as  the  Norwegian  fjords,  where  a  fresh  layer  often  overlies  a  more  saline 
layer.  Ships  caught  in  dead  water  were  held  to  the  internal  wave  speed 
unless  they  could  break  free  by  suddenly  increasing  their  speed.  Such  a 
sudden  increase  is  generally  possible  for  motor  ships  but  not  for  sailing 
ships. 

At  the  suggestions  of  Nansen  and  Bjerknes,  Ekman  [1906]  undertook  an 
extensive  study  of  dead  water.  Nansen  became  interested  in  the  subject  when, 
in  1893  off  the  Taimur  Peninsula,  the  Fvam  was  caught  by  dead  water  that 
reduced  the  ship’s  speed  from  a  normal  5  knots  to  1  knot.  Ekman  cited  this 
and  many  other  examples  of  dead  water  in  an  historical  summary,  and  he 
conducted  a  series  of  carefully  organized,  quantitative  experiments  which 
proved  conclusively  that  the  high  resistance  was  due  to  internal  waves . 

In  these  experiments,  ship  models  were  towed  in  tanks  containing  a  fresh 
water  layer  over  a  salty  layer,  and  the  towing  resistance  was  measured  as 
a  function  of  such  parameters  as  speed,  layer  densities  and  depths,  and 
hull  shape  and  size. 

The  atmosphere  is  also  stratified  and  internal  waves  frequently 
develop  in  the  lee  of  mountain  ranges.  Lee  waves,  or  mountain  waves,  are 
sometimes  made  visible  by  cloud  bands  at  the  wave  crests  which  remain 
stationary  with  respect  to  the  ground  despite  the  wind.  Reviews  of  the 
subject  have  been  given  by  Miles  [1969]  and  Turner  [1973] .  Although  a 
considerable  amount  of  theoretical  work  and  some  model  studies  have  been 
stimulated  by  the  atmospheric  problem,  the  application  of  the  results  to  the 
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present  problem  is  not  clear.  Model  studies  by  Davis  [1969],  for  example, 
do  not  show  good  agreement  with  theory.  The.  careful  model  studies  by  Ekman 
are  still,  despite  the  time  that  has  elapsed,  a  good  starting  point  for  an 
estimate  of  internal  wave  drag  on  pack  ice. 

Ekman  used  models  of  the  From  on  1/100  and  1/200  scales  as  well  as 
other  hull  shapes.  He  demonstrated  that  the  results  from  one  From  model 
could  be  applied  to  the  other  by  using  the  rule  of  Froude  for  ship  modeling. 
He  also  showed  that  the  results  could  be  scaled  to  give  results  in  close 
agreement  with  experience  for  the  full-scale  From.  Froude 's  rules  are 
based  on  dimensional  considerations  and  neglect  the  viscosity.  They  were 
originally  applied  to  surface  wave  resistance,  which  is  the  most  important 
drag  factor  at  higher  speeds  for  displacement  hulls. 

Ekman  extended  the  rules  to  cover  internal  wave  drag.  The  scaling  in 
the  case  of  wave  drag  is  in  terms  of  the  Froude  number  (F) ,  the  ratio  of 
ship  speed  to  the  speed  of  long  gravity  waves,  rather  than  the  Reynolds 
number  (R) .  The  drag  force  is  proportional  to  the  displacement  of  the  hull, 
or  cube  of  linear  increase  in  dimensions,  and  to  the  difference  in  specific 
gravity  between  the  layers.  This  allows  some  of  Ekman's  major  results  to 
be  summarized  by  a  plot  of  normalized  drag  force  versus  Froude  number,  with 
the  ratio  of  layer  to  keel  depth  as  a  parameter  (Fig.  1). 
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Fig.  1.  Normalized  internal  wave  drag  versus  Froude  number  with  the 
ratio  of  layer  to  keel  depth  as  a  parameter.  Curves  are  based  on 
Ekman's  [1906]  experiments. 
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Drag  results  for  homogeneous  water  are  also  shown  for  each  keel  depth. 
The  difference  between  the  homogeneous  water  curve  and  the  appropriate 
curve  in  stratified  water  represents  the  drag  due  to  internal  waves.  Two 
important  aspects  of  internal  wave  drag  are  shown  clearly  in  the  figure. 
First,  there  is  a  maximum  drag  when  the  Froude  number  is  0.73,  that  is,  when 
the  boat  speed  nearly  coincides  with  the  internal  wave  speed.  For  the  case 
of  keel  depth  less  than  layer  depth,  wave  drag  is  greatest  for  Froude  numbers 
less  than  one.  At  that  speed  the  largest  waves  are  generated  and  drag  is 
highest.  Internal  wave  drag  is  significant  only  between  Froude  numbers  of 
roughly  0.3  and  2.0.  Note  that  for  Froude  numbers  exceeding  0.73  there  is 
an  unstable  region  in  which  drag  decreases  as  speed  increases.  In  this 
speed  range,  ships  break  free  of  dead  water,  increasing  their  speed 
suddenly  until  limited  by  other  types  of  drag.  Second,  the  curves  of 
Figure  1  show  the  effect  of  changing  keel  depths  in  relation  to  the  depth 
of  the  fresh  water  layer.  Maximum  drag  occurs  when  the  keel  depth  is  the 
same  as  the  depth  of  the  fresh  water  layer. 

These  curves  may  be  applied  to  the  drag  of  pack  ice  provided  that 
layer  depth,  specific  gravity  difference,  and  displacement  of  the  ice  keel 
are  known.  Data  are  available  that  allow  one  to  make  reasonable  estimates 
of  these  parameters.  The  assumptions  are  that  the  shape  of  the  ice  keel 
resembles  a  boat  hull,  a  round-bottomed  one  in  the  case  of  Figure  1,  and 
that  the  Arctic  Ocean  stratification  can  be  represented  by  a  density 
discontinuity. 

Although  probably  only  the  keels  of  old  pressure  ridges  resemble 
the  hulls  of  ships  very  closely,  the  exact  form  may  not  be  particularly 
critical  for  an  estimate  of  internal  wave  drag.  The  different  hull  shapes 
that  Ekman  tested  showed  no  remarkable  differences  in  drag.  A  study  by 
Kovacs  et  al.  [1972]  of  an  old  pressure  ridge  with  a  keel  depth  of  13  m 
showed  a  smooth,  roughly  semicylindrical  underwater  form  somewhat  resembling 
a  ship’s  hull.  Younger  ridges  are  more  likely  to  be  composed  of  jumbled, 
unconsolidated  blocks  with  a  greater  form  drag  than  old  ridges,  although 
even  in  this  case  internal  wave  drag  may  not  be  affected  much. 

Although  stratification  in  the  Arctic  Ocean  is  actually  gradual,  its 
representation  by  a  discontinuity  is  a  reasonable  first  approximation  so  long 
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as  only  the  first  mode  is  being  considered.  A  density  profile  taken  at 
the  AIDJEX  main  camp  in  1972  (Fig.  2)  shows  a  low  salinity  mixed  layer  35  m 
deep  overlying  a  steep  density  gradient  extending  to  a  depth  of  300  m. 

Below  the  gradient  lies  nearly  homogeneous  water  extending  to  the  bottom, 
which  is  at  3800  m  in  this  area. 


Fig.  2.  Vertical  profile  of  temperature  (T) ,  salinity  (S) , 
density  (a^)  and  Vaisala  frequency  (N) ,  25  March  1972, 
72°07'N  149°00’W,  University  of  Washington  hydrographic 
stations  24  and  26. 
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Internal  wave  dispersion  curves  based  on  this  profile  are  shown  in 
Figure  3.  The  curves  are  based  on  a  matrix  method  for  calculating  internal 
wave  properties  in  a  rotating  ocean  with  many  layers  (Fliegel  and  Hunkins , 
in  preparation).  Note  that  over  the  nearly  level  portion  of  the  curve, 
phase  velocity  and  group  velocity  are  nearly  equal.  It  is  in  this  region  that 
internal  wave  drag  is  greatest.  The  appropriate  modeling  parameter  in  this 
case  is  the  phase  velocity,  which  is  about  200  cm/sec  over  the  level  part 
of  the  curve  for  the  first  mode. 


Fig.  3.  Internal  wave  dispersion  curves  for  the  first  four 
modes.  Calculations  are  based  on  an  18-layer  model  approxi¬ 
mating  the  density  curve  of  Fig.  1.  Solid  lines  are  phase 
velocity;  dashed  lines  are  group  velocity. 


A  two-layer  model  has  only  one  mode;  it  is  equivalent  to  the  first  mode  of 
a  multilayered  model.  For  the  two-layer  model,  the  interface  is  chosen  near  the 
center  of  the  gradient  at  a  density  of  1.026  g/cm3.  The  upper  layer  h  is  100  m 
thick  with  a  density  (p x )  of  1.024  g/cm3,  and  the  deep  layer  is  in  effect 
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infinitely  thick  with  a  density  (P2)  of  1.028  g/cm3.  The  phase  velocity  is 


C 


V 


1.028  -  1.024 

1.026 


x  983  cm/sec2  x  100  m 


195.8  cm/sec, 


where  g  is  the  acceleration  due  to  gravity.  This  is  sufficiently  close  to  the 
multilayered  model.  Some  experiments  that  Ekman  performed  on  multiple  layers 
showed  a  change  of  14%  at  most  when  the  density  change  occurred  in  several 
steps  rather  than  in  a  single  discontinuity. 

Information  on  the  numbers  and  average  depths  of  ridge  keels  comes 
from  U.  S.  Navy  nuclear  submarines  whose  upward-looking  sonars  give  profiles 
of  the  bottom  of  the  ice  [Lyon,  1961].  Actual  keel  depths  range  up  to 
47  m,  the  deepest  yet  observed  [Hibler  et  al.,  1972].  The  use  of  an  average 
keel  depth  appears  justified  for  rough  estimates  since  the  drag  is  approxi¬ 
mately  proportional  to  depth  over  the  range  of  interest.  For  the  central 
Arctic  Ocean  an  average  keel  depth  is  taken  to  be  10  m  with  a  spacing 
between  ridges  of  200  m.  The  ratio  of  layer  depth  to  keel  depth  is  thus 
10  for  the  Arctic  Ocean.  Since  the  largest  ratio  tested  by  Ekman  was  2.95, 
his  data  must  be  extrapolated. 

The  maximum  drag  appears  to  be  approximately  proportional  to  the 
reciprocal  of  this  ratio  for  the  range  of  interest.  For  a  layer/keel  depth 
ratio  of  1.14,  the  maximum  internal  wave  drag,  shown  by  the  double-headed 
arrow,  is  72  dyn/g  at  a  Froude  number  of  0.73.  This  is  multiplied  by  a 
factor  of  0.1  to  compensate  for  the  increased  ratio  of  layer-to-keel  depth 
in  the  Arctic  Ocean,  giving  7.2  dyn/g.  Assume  an  increase  of  200  in  linear 
dimensions  over  the  1/100  From  model:  the  actual  From  is  35  m  x  11.5  m 
x  4.4  m;  the  model  ridge  is  therefore  70  m  long,  23  m  wide,  and  8.8  m  deep 
and  travels  in  the  direction  of  the  long  axis.  The  70  m  long  ridges  would 
seem  to  occupy  a  large  fraction  of  the  pack.  Actually  a  shorter  hull  with 
the  same  draught  and  width  might  give  a  similar  internal  wave  drag.  In 
other  words,  the  internal  wave  drag  may  not  be  particularly  sensitive  to 
hull  length.  The  actual  From  displaces  800  tons,  and  so  the  model  ridge 
displaces  23,  or  8  times,  that  value.  Scaling  the  drag  we  have 
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D 

wave 


7.2  dyn/g  x  8  x  800  x  106  g  x  0.004 


1.84  x  IQ8  dyn 


Consider  an  ice  pack  with  one  of  these  ridges  spaced  at  200  m  intervals. 
The  internal  wave  stress  will  be 


1.84  x  ip8  dyn 
Twave  (200  m)2 


0.46  dyn/ cm2 


This  represents  an  extreme  upper  limit,  since  the  drift  speed  necessary 
to  attain  this  value  is  never  reached  in  nature.  This  maximum  value  of 
internal  wave  stress  is  achieved  at  an  ice  speed  of 

=  F  x  internal  wave  speed  =  F  x  /  g  x  sp.  gr.  diff.  x  layer  depth 

=  0.77  /  980  cm/sec2  x  0.004  x  100  m 


=  0.77  x  198  cm/sec  =  152  cm/sec  -  3  kt 

This  is  far  in  excess  of  actual  ice  speeds.  During  the  1972  AIDJEX  program, 

the  highest  observed  speed  was  26  cm/sec,  or  about  1/2  knot.  This  is 

equivalent  to  a  Froude  number  of  0.13.  Reference  to  Figure  1  shows  that 

internal  wave  effects  are  becoming  negligible  at  this  speed.  Ekman’s 

experiments  did  not  extend  to  such  small  Froude  numbers.  Extrapolation  of 

the  curve  suggests  that  the  internal  wave  drag  is  at  least  an  order  of 

magnitude  smaller  at  F  =  0.13  than  at  its  maximum  value.  Thus  under  actual 

conditions,  a  maximum  value  is  T  -  0.05  dyn/cm  . 

wave 

The  internal  wave  drag  must  be  compared  with  skin  and  form  drag  to 
determine  its  relative  significance.  In  terms  of  Prandtl's  boundary  layer 
theory,  the  roughness  parameter  value  at  the  AIDJEX  main  camp  in  1972  was 
zQ  =  4.17.  Assuming  a  reference  level  of  z  =  2  m,  corresponding  to  the 
depth  of  the  skin  friction  layer,  the  drag  coefficient  is  given  by 

CD  =  2k2/[ln(z/zQ)]2  =  0.021 

where  von  Karman's  constant,  k,  is  0.4.  The  stress  due  to  skin  friction 
during  the  period  of  fastest  drift  in  1972  was  then 

T..  =  pC„u2/2  =  1.026 x  0.021 x  (15)2/2  =  2.4  dyn/cm2 

skin  K  D  J 
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where  u  is  the  velocity  in  cm/sec.  There  will  also  be  a  contribution  from 
form  drag.  The  Reynolds  number  in  this  case  will  be 

R  =  p  ul/\i  =  1.026  x  26  x  103/10-2  =  2.7  x  106 

where  l  is  a  characteristic  length  and  y  is  the  viscosity.  Boundary  layer 
separation  has  occurred  for  most  geometric  forms  at  this  Reynolds  number,  and 
for  a  hemispheric  shape  the  drag  coefficient  will  be  in  the  neighborhood  of 
unity  or  perhaps  slightly  less.  The  form  drag  of  a  ridge  will  then  be 

D  =  p  £_  A  u2/ 2  =  1.026  x  1.0  x  23  m  x  8.8  m  x  (15) 2/ 2 

form  D 

=  2.33  x  108  dyn 

where  A  is  the  cross-sectional  area  of  the  ridge.  The  stress  will  be 

Tform  =  (2*33  X  1()8)/(200  ra>2  =  °*6  dyn/cm2 

These  calculations  indicate  that  internal  wave  drag  on  pack  ice  in  the 
Arctic  Ocean  is  negligible  in  comparison  with  other  types  of  drag.  This 
somewhat  justifies  the  neglect  of  internal  wave  drag  in  the  AIDJEX  field 
studies  where  only  skin  and  form  drag  have  been  considered. 

A  more  precise  estimate  of  internal  wave  drag  on  pack  ice  would 
require  direct  observations  in  the  field,  further  model  studies,  theoretical 
investigations,  or  some  combination  of  these.  Theoretical  investigations 
alone  will  probably  not  be  sufficient. 

The  study  of  internal  wave  drag  under  actual  conditions  in  pack  ice 
presents  difficulties.  The  internal  waves  that  produce  drag  will,  under 
steady  conditions,  have  crests  and  troughs  in  a  fixed  relation  to  the  ice. 

They  cannot  be  detected  from  a  single  observation  point  as  can  internal 
waves  generated  by  other  causes  which  travel  relative  to  the  ice.  Internal 
waves  moving  relative  to  the  ice  have  been  observed  from  drifting  ice  with 
thermistor  strings  and  with  echo  sounders.  However,  either  a  dense  horizontal 
net  of  fixed  sensors  or  a  mobile  sensor  is  required  to  observe  standing 
internal  waves.  The  unmanned  research  torpedo  developed  by  the  Applied 
Physics  Laboratory  of  the  University  of  Washington  would  be  an  ideal  mobile 
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sensor,  since  it  could  be  programmed  to  sample  along  paths  downstream  from 
ridges,  searching  for  the  undulations  in  temperature  and  salinity  indicative 
of  lee  waves . 

A  side-scanning  sonar  might  be  another  useful  technique.  It  has  been 
shown  that  acoustic  waves  with  a  frequency  of  100  kHz  are  reflected  by  the 
sharp  boundary  at  the  base  of  the  mixed  layer,  which  was  at  35  m  in  1972. 

A  side-scanning  sonar  with  high  frequency  might  be  able  to  detect  the 
undulations  of  this  interface  over  a  fairly  wide  area  and  give  a  three- 
dimensional  picture  of  it  in  time.  Further  model  testing  would  be  valuable 
to  extend  the  work  of  Ekman  to  specific  questions  about  pack  ice  drag.  Some 
points  to  be  clarified  would  be: 

1.  Drag  at  small  ratios  of  layer  depth  to  keel  depth. 

2.  Drag  at  small  Froude  numbers. 

3.  Effect  of  shapes  other  than  ships’  hulls. 

4.  Effect  of  continuous  stratification. 

The  calculations  have  been  based  on  two-layered  model  tests.  Although 
Ekman' s  results  seemed  to  show  otherwise,  there  is  still  a  possibility  that 
models  with  more  layers  will  show  high  drag  effects  for  large  ratios  of 
layer/keel  depth  at  certain  low  speeds.  In  a  two-layered  model,  only  one 
mode  of  motion  exists.  The  maximum  of  vertical  motion  is  at  the  interface, 
which  moves  up  and  down  as  horizontal  currents  in  the  upper  and  lower  layers 
flow  in  opposite  directions.  Higher  modes  with  more  than  one  maximum  of 
vertical  motion  in  depth  may  exist  in  fluids  with  many  layers.  The  higher 
modes  have  slower  phase  speeds  than  the  fundamental  mode. 

There  may  be  cases  in  which  a  keel  extending  a  small  distance  into  the 
mixed  layer  will  excite  higher  mode  internal  waves  of  large  amplitude  when 
it  moves  at  a  certain  lower  speed.  Ekman' s  limited  experiments  with  density 
gradients  could  have  missed  this  effect.  Note  that  the  fourth  mode  has  a 
phase  velocity  of  45  cm/sec  in  this  region  of  high  drag.  It  is  clear  that 
a  phase  speed  in  the  range  of  ice  speeds  will  occur  in  this  flat  region 
within  the  next  few  higher  modes .  These  modes  will  have  a  complex 
vertical  structure  which  might  be  excited  by  relatively  small  keels. 
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Vertical  current  shear  is  another  complication  which  could  be  of 
importance.  In  the  model  tests  the  fluid  has  a  uniform  motion  relative 
to  the  moving  object.  The  horizontal  velocity  in  the  ocean  does  vary  with 
depth  during  internal  wave  generation  by  pack  ice.  This  effect  has  been 
included  in  various  theoretical  studies  of  lee  waves  in  the  atmosphere. 

The  production  of  a  flow  with  horizontal  velocity  varying  in  depth  might 
require,  in  a  tank  model,  considerable  experimental  ingenuity.  Theory  may 
be  the  simplest  method  for  evaluating  the  importance  of  this  particular 
ef  feet . 

Although  there  are  still  a  number  of  unexamined  factors,  it  appears 
that  the  internal  wave  drag  generated  by  ridge  keels  is  not  likely  to  be 
an  important  factor  in  the  general  underwater  resistance  of  ice  floes  in 
the  Arctic  Ocean.  The  assumptions  and  extrapolations  involved  in  this  rough 
analysis  make  it  desirable  to  attempt  more  refined  calculations  and  experi¬ 
ments.  However,  the  neglect  of  internal  wave  drag  relative  to  skin  and 
form  drag  appears  at  present  to  be  justifiable  for  large-scale  ice  motion. 
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